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THE  TEMPERATURE  CALCULATION  OF  AN  ORTHOTROPIC  PLATE 
WITH  TRANSVERSE  SHEAR  TAKEN  INTO  CONSIDERATION 

S.  M.  Durgar’yan 
(Yerevan) 

1.  Statement  of  the  problem.  We  will  examine  a  rectangular 
orthotropic  plate  of  constant  thickness  h  referred  to  rectangular 
cartesian  coordinates.  We  will  combine  the  coordinate  plane  xoy  with 
the  average  plane  of  an  undeformed  plate,  we  will  direct  the  coordinate 
axes  along  the  principal  directions  of  elasticity  of  the  plate  mate¬ 
rial,  and  we  will  set  the  coordinate  origin  at  one  of  the  apexes  of  the 
plate.  The  deformation  factors  E^,  G^  and  the  thermal  coef¬ 

ficients  of  expansion  p^  will  be  considered  as  given  functions  of 
temperature  T,  changing  according  to  a  known  law. 

Here,  as  usual,  E^  and  p^  are,  respectively.  Young’s  modulus  art 
the  thermal  coefficient  of  expansion  along  the  coordinate  axis  lj  the 
Poisson  coefficient  characterizes  contraction  in  the  direction  j_ 
during  expansion  In  direction  JL;  the  shear  modulus  characterizes 
the  change  In  angle  between  principal  directions  i_  and 

We  will  assume  that  the  relative  deformation  throughout  the  plate 
Is  conditioned  by  free  thermal  expansion 
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-  3,T 


(1.1) 


Using  as  our  basis  a  theory  developed  by  S.  A.  Ambartsumyan  [1], 

we  will  consider  that  the  tangential  stresses  t„  and  r  can  be 

xz  yz 

represented  as 


T«  =  A  (0,  <P  (*.  !/).  T,(t  =  /.  (;)  l|>  (X,  y) 


(1.2) 


where  <p(x,y)  and  ^(x,y)  are  unknown  functions,  f^(z)  are  given  func¬ 
tions  which  characterize  the  law  of  change  of  tangential  stresses 
rxz  and  Txy.  throughout  the  plate;  fi(  +  h/2)  =  0. 

2.  Derivation  of  the  basic  equations.  Using  the  accepted  rep¬ 
resentation  of  tangential  stresses  (1.2),  having  integrated  over  z_ 
the  third  differential  equation  of  equilibrium  of  a  volume  element 
(assuming  absence  of  volumetric  forces)  and  having  designated 


we  will  obtain 


(2.1) 


y,  s)  =  «o(*.  y)  —  F\(t) Ft(t) (2.2) 

where  aQ(x,y)  is  the  integration  function  describing  the  distribution 
of  normal  stresses  along  the  average  plane  of  the  plate  (z  =  0)  and 
is  determined  from  the  conditions 

«,(*.  V,  ±h/2)~0  (2.5) 

Taking  into  account  (2.5),  from  (2.2)  we  obtain  the  value  of  the 
integration  function 

*  “• «  -  -r  ['■  (4)  +  '■  {-  4)]  £ + -T  ['•  (t) + '•  (-  4-)]  £ 


FTD-TT-63-529/1+2+4 


-2- 


and  also  the  differential  dependence 

which,  as  we  can  easily  see.  Is  the  third  equation  of  equilibrium  of 
the  plate  element  [1]. 

Using  the  found  value  of  the  Integration  function  aQ(x,y)  and 
having  set 

+  {i  =  1,2)  (2.5) 

for  stress  a  (x,y,z)  we  obtain  the  expression 
z 

«r.  *)  —  [«.— Fx (01  §+  i* -*<«)]§•  (2.6) 

We  know  from  earlier  articles  ([4],  pp.  74-75;  [1],  pp.  315-516) 
that  It  Is  advantageous  to  accept  the  functions  f±(z)  In  the  form 

— )  (2-7) 

Prom  (2.1),  by  virtue  of  (2.7),  we.  obtain 

f 

Thus,  taking  Into  account  (2.4) -(2.6),  we  will  have  c^  *  c^  -  crz 
(x,y,z)  -  0. 

It  Is  easily  ascertained  that  stress  o_  is  equal  to  zero  in  all 

z 

cases  when  the  functions  f^z)  and  fg(z)  are  even  functions  of  the 
coordinate  z  and  are  equal  to  each  other. 

From  this,  in  subsequent  calculations  we  will  set  a_  ■  0. 

On  the  basis  of  Franz  Neuman* s  hypothesis  [2],  from  the  gener¬ 
alized  Hooke  law  we  have 

s,  =  #u*x+  +  QiT,  ov  *=  Bi$ex  +  Dney  +  Q»T 

'  -Cxy  =  S"'XV.  V  “  D**V,'  XXt  “  D"?XI 
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where,  as  usual  (for  example,  [1]), 


Du  ~ 


1  “  l*xv*ly*  * 

E. 


Du  i 


1  -*WS/x’ 


On  «=  T-r 


-tryx 


1-^w*  1  -  ^UX 


— n  n  (3«  +  >‘»A)'  G»-  i  — 


(2.9) 


rafyryjt 

Du  -  Gv,, 


'  Gtx,  Du  =  G, 


*v 


Using  the  known  relations  between  displacements  and  deformations 

and  also  taking  into  account  (l.l),  (1.2)  and  the  last  two  equations 

of  (2.8),  for  projections  of  displacement  u  ,  u,  and  u  we  obtain 
v  x  y  z 


du> 


'*  “  “  ~  *  7*  +  ,Pr»*  uy"  *  —  *  ^g*  —  ly  +  Wi,  utm.u+r 


(2.10) 


Here 


M  “  M*  l*«0'  9  “  ttv  ll*0»  W  “  Ml  !*— 0 

*  *  , 
.  r,~.  ,  r*t 

’•  v*.  y.  *}  *■  44  —  v*.  v.  *j  *■  \ 


ra.u) 


;  •%(«*  . 
1  '*« 


.  c;--ai..AS 


stress-1 


.  <r. 


\ 

) 


ne  .••-r.v*r  *  -n-jctulaa  it  .Is  not  di:."- 

-j.r.0  ,viw*  :•  means  of  the  ra*JU 


unknowns  u,  v,  w,  9,  and  $ 

Tx~  f  v,*  =  c“i7+c“i7“JC,‘S-”x‘,lF+Cm*9  + 


+  Cm^  +  Cn,  ^  4  f  m  -  *»“  -  *‘*V  +  ®» 


”1V,“ 


C^  +  C“S-X”S-^^+C“‘9  + 


+  Ciav'f  4  Cm  4  Cm  —  Ri*>x  4  ®» 

Smm  \  T*vdl  “  C"  (|^  4  4r)  ~ ^ S^Sy  +  c,,iV,f  + 
-Vi 

4  C*m*^  4  Cui  -^-  4  cw^  —  ^vo*  “  "*>V 

T  ^-^(-r)*  f  v'*-^(4-) 

-v«  -i« 


4  *mv*  4  ^ui  "gj  4  ^  —  Bin*  —  Bin*  4  811 


Af„-  1j*  «„**- 


*»  ir + ~  ^  ^  ,5’ + *m*9 + 


4  *mv*  4  *m  H 4  JC»  |j£  -  V  -  *mv 4  Ou 

-i/» 

4  4  4  "aj  ^*i*  "xiy 


Here  we  use  the  designations 


(-  T  V*.  Jf,  —  T  ct-  T  C„-  ^  a,iyr 

— ' ifl  —A/I  —Art  —Vi 

Y  <y_  V  »,45t*.  V-  f  ».-£- 

— 1 Art  —a/*  —Art 

(<-11,  12,  22,  66;  /-l.  2;  v  =  »,  y) 

^  (•'.*»  *. v;  /-y.  *;  v-o,  i;  «+/) 

—i/I 

«’•  T  8i,W'il  (i  —  11,  12,  22;  v  -  y;  /  -  0,  1) 

,  y  \  *  jv 


(2.13) 


(2.14) 


Art 

8«-  (  Qt T*1*  (*  =1,2;  /  —  0, 1) 
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Introducing  (2. 13)  into  the  equations  of  equilibrium  of  the  plate 
element,  for  determining  the  five  unknowns  u,  v,  w,  9,  and  f  we  obtain 
the  following  system  of  five  differential  equations  with  variable  coef- 
ficientsi 


For  convenience,  the  designation  ^--37- -®  is  used  in  (2.15)- 
(2.19);  from  this. 


J.  Certain  particular  cases.  The  case  when  the  temperature 
changes  only  throughout  the  plate,  T  =  T(z) .  In  the  examined  case  the 
system  of  equations  (2.15)-(2.19)  is  simplified  and  reduces  to  the 
following  system  of  five  uniform  differential  equations  with  constant 
coefficients 


(Cu  •gj  +  Cu  *§7)  “  +  (Cl,  +  Cm)  gjgy  —  [jfu  +  (Z«  +  2X**)  dxd\j*\  w  + 

1  a*  a*  \  a** 

+  (Cm  gxT  +  Cmi  -57]  <P  +  (Cm  +  Cm)  giity  “  0 

(C,t  +  Cm)  -^  +  (Cm  ^•  +  Cn-^i)»- [(*«  +  2A'«)  aSy +  "  + 

+  (Cm  +  Cmi)  jjfij  +  (Cm,-|?  +  C«  ♦  -  0 


do  aik 

1T+ ajr 


o 


(30) 


(*u  I?  +  *m  -J7  )“  +  (*«  +  +  (fl»  +  2fl«)  §^i]  »  + 

[  a*  a*  /  fc  \T  •  a*'i> 

*rn  a?  +  ^Mi  -57  -  2^  (x)J  *  +  aiaiT  ”  0 

(K„  +  K")  +  (^m  a*T  +  *«,  37)  »  —  [(Ait  +  2Am)  foZfoj  +  Am  ^7]  w  + 

+  (*m  +  K*u)  +  [A'm  a7  +  Ktn  -57—  2^  (”)].'♦  “  0 


The  case  when  the  temperature  changes  only  along  the  coordinate 
axes  x  and  y,  T  -  T(x,y) .  In  this  case  we  have 


Bi-Bi(*.  v),  Qj*=Q>(*.  y) 

/•-P.t*.  /,  =  i 


(i  =  11,  12,  22,  44,  55,  06;  ;  -  1,  2) 


£(gj) 

dx 


'«  2  *  ay 


(3.6) 
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Using  (2.7) »  from  (2.1)  and  the  last  two  equations  of  (2.1J.)  we 
obtain 


„  1  //>*  **\  „  1  I  h*  **\  „  /  /.  \  h» 

vt-tv*  (3.7) 

Taking  into  account  (3.6)  and  (3.7) >  from  (2.14)  we  will  have 

b{  A* 


h* 


^“p^rss*  V“B‘  77(1^)125 

(*'  —  11.  12, .22.  66;  v  -  y) 

„  ,  _  aa  (3,T)  r¥  d=>(0,T) 

*•  “  24  dv*  '  "»  “  2?  Jvd/-  •  *U 

(i  —  11,  12,  22;  V  —  y;  /'  *»  y,  *) 

8le-QiTA.  6tl-.0  (»-1.2) 


(3.8) 


Inserting  (3.8)  into  (2.15) -(2. 19)  we  obtain  the  following  system 
of  differential  equations  with  variable  coefficients 


B„h*  a*(P,T)  (pI,  +  2BM)At*’(3rT> 

~d*^  +  54 - a^r 


24 


g(QiT) 

09 


dy)  +  dy  (*>*  a*)]“  +  [^(B“^)+^(Z,»^)]*- 


J}nh*  ^(3,T)  (P„  +  2D„)  /.»  y(3,T) 

dy*  ■**  24  dxldy 


24 


d<p  d\b 

~oi  +  V  =  0 


a(Q«T) 


+ 

+ 


~[a*  (c“  ^  +  c>«^r)  +  2-^(fl«  a5^)] 
tfi>  £  -k  {ik  ??)] + T5-  iy  [B-  (i  Jp)] 


»  + 


10  {  a*  [B“  dy  ( p44  ay- )]  +  i  [°M  37  (ifcr  3F-)]}  ♦  -  0 


<P  + 


r, 5 

a3  \ 

_a_. 

( 

a* 

“L2  a* 

[o„ 

dxdy  ) 

dy 

[Bn 

a** 

+  Bn 

r  3  r 

_d_ 

/  1  a*  \ 

iT 

a 

r 

a  / 

+  75' 

la^r'4 

dy 

[liu  dy*  1 

0+ 

dy 

(*» 

ir(: 

/**  a  1 

‘  d  1 

’  1 

a*  \1 

h* 

a  1 

r 

1 10  '  a*  [ 

*"1 T\ 

1 

a?)J  + 

10 

dy  1 

[Bn 

aJT  (' 

u>  + 


■  + 


(3.9) 


(3.10) 


The  system  of  equations  (2.15) -(2. 19)  was  broken  down  into  two 
systems:  (3.9)  (plane  problem)  and  (3. 10)  (bending  without  expansion- 
contraction  of  the  average  plane) . 
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4.  Cylindrical  bending  of  a  plate.  Having  directed  axis  x  along 


the  span  of  the  plate,  from  (2.15)-(2.19) ,  taking  Into  account  (2.7), 
we  obtain 


<Pw  \  d  ( 
dx *  7  "t  dx  V 


C,u  lx  +  Cl 


9 


£  (W — Qi») 
dx 


d<f/dx  =  0 


d_ 

dx 


lx  +K"* 


d  (Jim* — 8u) 

dx 


(4.1) 


From  the  second  equation  of  (4.1)  we  find  ?  *  const.  Considering 
this,  and  having  integrated  once  the  first  and  third  equations  of 
(4.1),  we  will  have 


du 


dSo 


C\\  -r-  —  Km  =  —  Cmx<p  +  Riux  —  6io  +  Ci 


(4.2) 


(4.3) 


where  C^  and  C g  are  constants  of  integration. 

For  definiteness,  having  taken  C^D^  -  K^2  ?  0,  (4.2)  and  (4.5) 
we  will  obtain 

■  CiiDu-Kn'  [*  (*»*»»*  -AlCm*  -  A'u  -55-)  + 

4*  C\Du  —  C%K\\  +  DiiRm»x  —  K\\Rui *  —  Au0 10  4" 

*£  1  L/c  JC  x  KCx  c  —)4-  (4.4) 

°  c uii  11 — a’u*  L9  -c“  u  ;  + 

4-  CjXn  —  C,Cu  4-  KiiRuo*  —  Culim*  —  An0i«  4-  Cu0iij 

Integration  (4.4)  does  not  present  any  great  difficulties  for  the 
known  law  of  temperature  change  and  for  given  functions  which  express 
the  dependence  of  the  deformation  f&otor  and  the  thermal  coefficient 
of  expansion  on  temperature. 

The  six  constants  of  Integration  which  enter  into  expressions  9, 
u  and  w  are  determined  from  the  conditions  that  the  plate  is  attached 
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along  the  sides  x  -  0  and  x  ■  l.  [There  are  three  constants  of  integra- 
on  <p,  C1  and  C2  in  (4.4)  and  three  also  appear  as  a  result  of 
integration  of  (4,£).] 

Example  of  ths  oaloulatlcn.  We  will  examine  a  otae  when  aero 
temperature  is  maintained  at  the  lower  surface  of  the  plate,  while 
at  the  upper  surface  the  temperature  changes  linearly  from  0°  to  400°, 


i.e.. 


T-  400  y  «fa«n  i-.-2.T-0  te>— j- 


(4.5) 


We  can  easily  see  that  the  temperature  function 


t  -  «»t  (t  +  t)  ^ 

satisfies  the  known  equation  of  heat  transfer  and  the  boundary  con¬ 
ditions  (4.5) . 

We  will  take  [5] 

&  c 

—  -j?  —  1  —  0.0005T,  p,  -  p,  -  p  -  eowt, 

Ms*  —  |»i  —  const, '  |i^  —  m  —  oonat  ' 

Introducing  (4.6)  into  (4.7)  and  introducing  new  variables 


2*  ’  x 

“ — nsr 


(4.8) 


from  (2.9),  (2.14)  and  (4.4)  we  will  obtain  the  values  of  the  first 

derivative  of  displacement  u  and  the  second  derivative  of  displace-  • 

ment  w  with  respect  to  variable  <0,  which  contain  logarithmic  terms. 

Takinf?  Into  account  the  boundedness  of  variable  co 

♦ 

0<«<i/i0 

i  -•  ■  *  ,'ij  under  the  condition  l  £  x  £  0)  and  given  the  accuracy  of 

"l..:.  ,  we  can  substantially  simplify  the  obtained  results. 

■.n-  '.-'v?  expanded  the  logarithmic  functions  into  power  series  and 


- •  .  j.  —.j?.  j  1 1+ 2+4 
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It 

having  limited  ourselves  to  an  accuracy  of  1  +  o>  »  1,  after  integra¬ 
tion  we  will  obtain 

u  —  u  (oo,  <p,  Ct,  Ct,  Ct),  (»,  9,  Cl,  C»,  Cu  Ct)  (4.9) 

Introducing  (4.2)  and  (4.9)  into  (1.2),  (2.10),  (2.12),  and  (2.13)* 
we  can  easily  find  the  stresses,  displacements,  forces,  and  moments 
expressed  by  constants  q>  and  (i  ■  1,  2, . . . .,  5) ,  the  values  of 
which  are  determined  from  the  conditions  that  the  plate  is  attached 
at  the  edges  x  =  0  and  x  ■  l  (correspondingly,  cd  -  0  and  cd  -  0.1)  . 

We  will  examine  a  case  when  along  the  sides  x  «  0  and  x  «  l  where 
«  z  «  o,  the'  following  conditions  are  satisfied. 


®u. 


(4.10) 

From  (2.10),  (4.9),  and  (4.10),  taking  into  account  the  bounded¬ 
ness  n(  |ti|  £  1),  we  will  have 


0, 

92.59  + 0.1481ft* 
C, - T—f - 


3S7.4  S* 

c,— 1 


(4.11) 


Here 

b  ft*  ..  „P(i  +  mi 
h  ~  TT-  ltr££r 

Obviously  from  (4.11),  in  the  case  h*  -  0.114-1.00,  the  correc¬ 
tion  resulting  from  taking  into  account  transverse  shear  when  deter¬ 
mining  Cg  can  reach  5-44jf. 

Introducing  (4.4)  into  (2.13)*  for  the  deflecting  moment  Mg  we 
will  have 


Mt  m  0.06339ft**  +  Ct 
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from  which  it  is  evident  that  talcing  into  account  the  influence  of 
transverse  shear  can  in  certain  cases  very  substantially  influence 
the  calculated  values  (for  example ,  the  value  of  the  deflecting  moment 
In  oress  Motion  x  -  0  in  the  mentioned  example) . 

Submitted  July  10,  1961 
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DISTRIBUTION  LIST 


DEPARTMENT  OF  DEFENSE 

Nr.  Copies 

MAJOR  AIR  COMMANDS 

Nr.  Coplea 

AFSC 

SCFDD 

1 

DDC 

25 

TDBTL 

5 

HEADQUARTERS  USAX 

TDBDP  . 

TDENP  (Roney) 

2 

1 

AFCIN-3D2 

ARL  (ARB) 

1 

1 

AEDC  (AIT) 

AFITC  (ITT) 

ASD  (A3YIM) 

1 

1 

2 

BSD  (BSF) 

1 

BD  (1ST) 

1 

8SD  (SSP) 

2 

OTHER  AGENCIES 

CIA 

1 

NSA 

6 

DIA 

9 

AID 

2 

OTS 

2 

ABC 

2 

PUS 

1 

NASA 

1 

* 

ARMY  (FSTC) 

3 

NAVI 

3 

NAFBC 

1 

POE 

12 

RAND 

1 
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